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Abstract 

We consider the fermionic bound states associated with a soliton-antisoliton pair in 
1+1 dimensions which have zero energy when the solitons are infinitely far apart. We 
calculate the energies of these states when the solitons are separated by a finite distance. 
The energies decay exponentially with the distance between the soliton and antisoliton. 
When the fermion mass is much larger than the boson mass, the energy simplifies sub- 
stantially. These energies may be interpreted as a contribution to the effective potential 
between the soliton and antisoliton. The character of this contribution depends upon 
which fermionic states are occupied. Performing the analogous calculation for the sim- 
plest (3+l)-dimensional soliton system, we find no fermionic energy shift. 
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The (l+l)-dimensional Dirac equation, when coupled to a scalar field, has the in- 
teresting property that with the scalar field in a single-soliton state, the fermionic spec- 
trum can include a single zero-energy mode, so that the electric charge of the system 
is half-integral [1]. However, while the one-soliton case is exactly solvable, the problem 
of fermions in the presence of multiple interacting solitons is difficult to analyze exactly. 
Graham and Jaffe [2] have developed a numerical method for calculating the fermionic 
energies in a multi-soliton background. We shall examine the two-soliton case using an 
alternative method, deriving an analytic approximation for the fermionic energies. Our 
approximation method can also be applied to problems in 3+1 dimensions, where numer- 
ical calculations are difficult. 

The free Dirac equation in 1+1 dimensions is i^d^ip + rmfj = 0, where a suitable 
choice for the Dirac matrices is a = a 2 and (3 = o"i. If we use this representation, then 
the Hamiltonian H = ap + (3m becomes 



H = —ia 2 d x + a±m 



— d x + m 
d x + m 



(1) 



We shall replace the fermion mass m with a scalar potential U (x) , which is time- indepen- 
dent (or varies very slowly in time), so that energy is conserved. 

We are interested in the case where the potential U(x) that appears in the Dirac 
equation is derived from a scalar field <j>. We shall take the Lagrange density for this field 
to be 

£ = l<j> 2 -l(d x <t>f-V(4>). (2) 

We shall concentrate on the "Mexican hat" 4 potential, 

although similar results hold for for the sine-Gordon potential [3, 4]. 

The 4 potential has two degenerate vacua, at = ±f . There are also solutions that 
interpolate between the two vacua (i.e. solitons). For a stationary soliton centered at 
X — Xq } the solution is 



^s,x (x) = -tanh 



(4) 



where /i 2 = 2a 2 is the mass squared of the light quanta of the scalar field. As x — > ±oo, 
4>s,x { x ) —> if- There is also an antisoliton solution 4>a,x { x ) — ~4 > s,x (x), which goes 
the other direction, from +f at x — — oo to — | at +oo. Finally, there are solutions that 
resemble spatially alternating solitons and antisolitons. However, treating these solutions 
as collections of individual solitons and antisolitons is a good approximation only when 
the solitons and antisolitons are widely separated, so that they disturb one-another only 
minimally. In general, these solitons must be moving relative to one-another, but we shall 
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assume that this motion is very slow, so that <p{x) depends very weakly on time. In three 
spatial dimensions, there do exist exact multi-soliton (multi-monopole) static solutions to 
bosonic field equations. 

We may introduce massless fermions coupled to this scalar field by a Yukawa interac- 
tion, U(x) = g<p{x). When is in a vacuum state, this generates a fermion mass m = g^. 
When cj) has the soliton profile <fis,x , there is known to be a zero-energy bound state 
localized around the zero of the soliton [1, 4, 5], with wavefunction 



ip + — C | cosh 



^(X-X ) 



-2m/ )x 



1 





(5) 



C is a real normalization constant. For the antisoliton 
wavefunction is 



1>- 



C { cosh 



— \x 



x 



A, xo the corresponding zero mode 

2m/ /i 





1 



(6) 



We shall investigate how the energies of the fermion bound states are affected in the 
two-soliton case. We shall consider an antisoliton centered at x\ and a soliton centered 
at X2, where /i(x 2 — xi) 1- We shall approximate as a sum = (f)A, Xl + 4>s,x 2 + f an d 
use vp + and ip- to represent the bound state wavefunctions (5) and (6) centered at X2 and 
xi, respectively. 

When the solitons are far apart, we may expect the potential due to one soliton to be 
a small perturbation of the Dirac equation for the bound state localized around the other 
soliton. In this spirit, we shall calculate the energy expectation values for the symmetrized 
states ^± = ^75 ("0+ ± ip-)- For Hip + , we have 



Hip+ = [-ia 2 d x + a 1 g<f> StX2 ]ip+ + o x g 
= + a 1 [g(f) A ,x 1 +m\ip + . 



a 



(7) 



We shall first evaluate the energy expectation values E± = /f^ dx W±H^± under the 
approximations ^ 3> 1 and yu(x 2 — x\) 3> 1, where elementary expressions can be obtained 
explicitly. Then we shall present an exact analysis for general ^, which reduces to the 
previous when 3> 1. 



The characteristic 



To find E±, we need to integrate quantities such as (tp + aiijj-)(f)A,xi 
size of 4>a,xi is ~> an d the characteristic decay length of ip± is ^. If ^ ^> 1, then the 
fermion wavefunctions decay very little over the width of the antisoliton, and it is a 
good approximation to replace <Pa,xi by the step function — |sgn(z — x\) in the integral. 
Approximating <f>A,xi with this step function gives 



Htf) + = <7im[l — sgn(x — xi)]ifj. 



The corresponding expression for Hif)^ is 

= <7im[l + sgn(x — a^)]^-- 



(8) 



(9) 



2 



Combining (8) and (9), we get 



HV± = -^=<7ito{[1 - sgn(x - xi)]if; + ± [1 + sgn(x - x 2 )]ip-}. (10) 
v2 

Since ip\aiip + = ^_a\ijj- = 0, only the cross terms in ty±H^± are nonzero. Using the 
relations 



ip+<7iip- = V ; - cr i'0+ = C 2 \ cosh 
the expression for E± is 



(X — X\ , 



cosh 



x - x 2 ) 



-2m/ fi 



;n) 



E± = ±m 



dx < cosh 



x 



2 — sgn ( x — y 



2 

Ax\ 



2 / 



Ax\ 
I + sgn I x — y — 



^[x-y+—j 



cosh 
Ax 



Ax\ 



!(-'-¥)]} 



-2m//j, 



;i2) 



where y = a2 + Xl , and Ax = x 2 — xi is the separation between the solitons. Shifting he 
integration x — > (x — y) gives 



i£± = ±m— / dx { cosh 
2 J —00 



x + 



Ax 



x 



Ax\ 

2 — sgn [ x H — — I + sgn [ x — 



cosh 
Ax' 



Ax 



x 



-2m/ n 



(13) 



which depends only on Ax. Thus we see that the zero-energy eigenstate of the one-soliton 
background bifurcates into two states, symmetrically displaced above and below E = 0. 
To proceed further with our estimate of E±, we must approximate the wavef unctions, 
using 



cosh 



-(x-x ) 



6{x - x )e^ (x " Xo) + 6(x - x)e-^ x ~ Xo) 



(14) 



The approximation (14) is natural in this context, since the approximate wavefunctions 
produced by (14) are the exact wavefunctions corresponding to perfectly rectangular soli- 
tons. Using (14), E± reduces to an integral over exponentials and step functions, which 
may be evaluated exactly, yielding, 



E± « ±(16) m/At CV 



■mAx 



(15) 



It only remains to calculate C. Using the same approximation (14) for cosh and demanding 
that / ip+ip+ dx = 1 gives 



C 2 « (16) 



-m/fj, 



so that 



E± « ±me 



m, 



-mAx 



(16) 
(17) 
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An interaction energy of this nature has been previously suggested [1]. 

We are now in a position to show that the ty± are the correct approximate states 
for which to calculate the energy. If we instead used some other linear combination 
{bip + + c0_), with \b\ 2 + |c| 2 = 1, the expectation value of the energy would be E' = 
2(Re{bc*})E + . E' takes its minimum value (corresponding to the best approximation to 
the lowest-energy wavefunction) when c = —b. The best approximation to the higher- 
energy wavefunction must be orthogonal to our expression for the lower-lying state, and 
so has c = b. So the best linear combinations of -0+ and vp- are 4^ (vf) + ± V 7 -)- 

To obtain the result (17), we needed to make the approximation ^- 1. We now 
present an exact evaluation, which shows that the above is the leading term for large — . 

If we insert the exact expressions (4), (5), and (6) into (7) (and the corresponding 
expression for Hi^J), we get 



± 



Gllll 



tanh 



V2 

±—=ovm < 1 + tanh 
y/2 I 



x - y + 



Ax 



Ax\ 
~2~) 



(18) 



Multiplying this by ty± from the left and integrating over the shifted variable (x — y) gives 

r 

J—c 



E. 



±m- 



C 2 



dx \ sech 



x 



tanh 



2 ' ■'" 



a ( Ax 

2\ Z+ - 
Ax 



sech 



+ tanh 



2 ' " 



ji ( Ax 

Ax 
~2~ 



2m/// 



(19) 



When the integral is performed, the two tanh terms contribute equally, giving 

2m /fi 



E. 



± 



±mC 2 
x < 1 + tanh 



/ dxlsech 

J — CO I 



■fi 


( Ax\ 


sech 


■ft 






(-£)] 


2 






2 





fi ( Ax 

— ix 

2 V 2 



We may simplify (20) using hyperbolic identities, in particular 

sech a sech b 



sech (a + b) 



1 + tanh a tanh b 



(20) 



(21) 



Letting a = |(x - 4f ) and b = f Ax in (21) gives 



E± = ±mC 2 sech™ ( —Ax 
V n 



sech 2 ™ 


m 
n 




' {l+tanh[= 




]} 


|l + atanh 


[?(*"¥)]}" 



, (22) 
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where we have introduced the dimensionless parameters n = — and a = tanh(yAx). 

Making the substitution u = tanh[^(x — ^)] in (22) simplifies this integral substantially, 
giving 



E ± = ±mC 2 sech n Ax j J ^ 



' l , ra (1 + u) n (l -u) 
du 



n-l 



(23) 



m (1 + cra) n 

To get a final expression for the energy, we must again calculate C. The normalization 
gives 

1 



C 2 



dx sech 



2n 



in 



(x - x ) 



= r 1 du-n-u 2 ) 

J-i m 

n j— F(n) 
V 7r i 



n 

2\n-l 



m T(n+|) 



(24) 



which agrees with (16) as n — > 0. 

The remaining integral in (23) may also be found analytically. The result involves the 
hypergeometric function F(ai, a,2',bi;z). In terms of this F, the integral is 



1 , (1 + u)"(1-m) 
du - 



n-l 



(1 + au) n 



1 (n + oj v 1 + ct 



(25) 



Fortunately, this simplifies significantly when combined with (24), so that the final formula 
for the fermion energies is 



E± = ±m sech™ (—Ax) (1 + a)~ n F (n, n; 1 + 2n; . 

\n J V ' 1 + aJ 



(26) 



Since the solitons are far apart, a ~ 1 — 2e 2 ™ A:E = 1 — e is close to one, and it is 
natural to expand (25) around a — 1. If we expand the integral in (25) to 0(e), we find, 



1 , + 
du — 



n-l 



1 + aw)™ 



J ^duil-u) 11 - 1 + e n J 



1 , u(l-u) 
aw 



n— 1 



1 + tl 



(27) 



The 0(e°) term in (27) has value However, the (9(e) term is divergent, because 
F(n, n; 1 + 2n; z) is not analytic at z — 1. Instead, we have 



1 , (1 + m)"(1-u) 
cm - 



n-l 



2™ 



+ n2™" 1 (elne) + C(e) 



(1 + au) n 

Using (28) and e = 2e~^~ Ax , the O(elne) expression for E± is 



(28) 



2™ T(n + |) 
7^T(n + l) 



sech™ ( — Ax 
n 



1 — n (2mAz) e 



(29) 
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As n — > 0, this agrees with the earlier result (17), through terms of G(n). 

The integral (25) may be evaluated in terms of elementary functions whenever n is an 
integer. Of particular interest is the case n — 2, which corresponds to a supersymmetric 
Lagrangian. For that case, the integral is 



/. 



. (1 +«)'(! 

(1 + au) 2 a 4 



-4a 2 + 6a) + (a 2 + 2a - 3) In 



a 



1 — a 



(30) 



and the energy is 
3 



E. 



± 



±-msech 2 
2 



rn 



Ax 



a" 



{-2a 2 + 3a) + (^ Ax ) (« 2 + 2a - 3) | 



(31) 



The first term on the right-hand-side of (30) gives the O(e ) term, while the second term 
is C(elne). 

The expression (26) gives the energies for the fermion bound states that lie closest to 
zero energy. When n < 1, the one-soliton system has only the single bound state (5) [5], 
so the soliton-antisoliton system we are considering has two bound states, with energies 
given by (26). The fermion energies depend on Ax, so (26) generates a contribution to the 
effective potential between the soliton and the antisoliton. The character of this potential 
term depends on which fermionic states are occupied. In the vacuum, the state ^_ is 
occupied, while ^ + is empty. E_ is negative and for a 1 becomes more negative as Ax 
decreases, creating an attractive interaction between the soliton-antisoliton pair. Since 
this interaction exists in the absence of any fermions or antifermions, it must be generated 
by virtual particle effects. If either a single fermion or a single antifermion is present (i.e. 
both states are either filled or empty), the total energy vanishes, and these states do not 
contribute to the effective potential. Evidently, the presence of a single particle causes 
a cancellation of the virtual particle corrections mentioned above. In the presence of 
both a fermion and an antifermion occupied and empty), the potential term 
is repulsive. This presents an interesting picture, with the fermions introducing novel 
interactions between the soliton-antisoliton pair. 

If the positions of the soliton and antisoliton are reversed, so that x\ > x 2 , a calculation 
of the energies proceeds along the same lines. We could extend (26) by replacing Ax by 
|Ax| (both in the sech term and in a) and multiplying by sgn(x 2 — Xi). The energies 
obtained would be valid whenever /x|Ax| ^> 1. However, there are additional subtleties 
that arise when we allow Ax to change signs. This involves the solitons passing through 
one-another. The 4 solitons do not pass through one-another to emerge undistorted [6] ; 
a sine-Gordon soliton and antisoliton will pass through each other without distortion, 
but when the solitons are in precisely the same place, the field takes a vacuum value 
everywhere, so there are no bound states at that instant [3]. In either case, it is impossible 
to identify which bound state for x 2 > x\ corresponds to which bound state for x\ > x 2 . 
However, (26) does have the attractive feature that the energies E± approach the energies 
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of stationary continuum solutions as Ax — > (even though (26) is not valid for small 
values of the separation). 

The simplest (3+l)-dimensional soliton with a known zero-energy fermion mode is 
the 't Hooft-Polyakov monopole [1, 7] coupled to isospinor fermions (in the fundamental 
representation). The monopole is described by a scalar field $m = an d a vector field 

A° M = 0, A l M = A l ^T a , where T a is the 2x2 isospin generator. Using the representation 



(T 

a 



-/ 
/ 



(32) 



for the Dirac matrices, the associated zero mode is 



1>. 



M 



x+ 





X? 



/(|r-r M |)( S ;^- S ^), 



(33) 



where v is the spinor index, n is the isospin index, ym is the position of the monopole, and 
/ is a known function of the boson fields. The s± n are orthonormal spinor, isospinor basis 
vectors. The corresponding zero mode for an antisoliton with boson fields 3% A° A = 0, 
and A\ is 





X- 



X- 



f(\r-r A \)(s» + s n _-s»_sl). 



(34) 



To apply our method to this system, we consider the approximate monopole-antimono- 
pole configuration $ = ^m+^a+^o (where $o is the vacuum value of the field as r — > oo), 
A^ = A^ M + A\ and calculate the expectation value of the energy for the state 



± 



1 

V2 



X+ 

X- 



(35) 



in this background. This energy is 



E 4 



d s r [xl X ] 



cr-p- a--A a T a -iG$ a T a 

o--p- (T-A a T a + iG$ a T a 



X+ 
X- 



(36) 

where G is the strength of the Yukawa coupling. Since the states x± are singlets of spin 
plus isospin, the terms with p and $ vanish immediately. The vector field term gives a 
contribution proportional to 8fA ta , which vanishes for the monopole-antimonopole profile, 
because A m is traceless in (i,a), since it involves e MJ . So in this approximation, there are 
two degenerate zero-energy fermion modes for the monopole-antimonopole system. 

The isovector fermion does not obviously produce zero-energy modes. Its properties 
are under investigation. 
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